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INTRODUCTION 


A  quantitative  knowledge  of  solid  propellant  response  functions,  both 

pressure-coupled  and  veloci  ty- coupled,  is  a  critical  input  parameter-  for 

rocket  motor  stability  analyses.  There  have  been  several  experimental 

techniques  developed  to  measure  propellant  pressure-coupled  response  functions 

(Refs.  1  &  2).  The  most  widely  accepted  method  for  propellant  response 

measurement  i.3  the  T-burner  method.  The  T-burner  has  solid  propellant  placed 

at  the  opposing  endface3  of  a  center— vented  tube.  Pre33ure  oscillations 

develop  between  the  two  propellant  3urface3  at  the  fundamental  longitudinal 

frequency  of  the  tube.  Pressure-coupled  response  data  is  derived  by  comparing 

the  measured  rate  of  growth  and  decay  of  the  pressure  oscillations.  The 

T-burner  i3  limited  to  one  pressure  oscillation  frequency  for  each  chamber- 

geometry.  A  complete  investigation  over  a  large  range  of  frequencies  requires 

many  different  tube  lengths.  Other  methods  of  measuring  the  propellant 

pre3sure-coupled  response  incorporate  externally  driven  pressure  oscillations 

in  order  to  examine  an  entire  frequency  range  with  a  3ingle  burner  geometry. 

One  3uch  method  utilizes  the  rotating  valve  burner .  Thi3  burner  produces 

pressure  oscillations  by  using  a  rotating  valve  to  periodically  modulate  the 
0 

exit  area  of  the  exhaust  nozzle.  The  pressure  oscillations  within  the  chamber 
are  a  function  of  the  propellant  response  function  and  the  gas  dynamics  of  the 
combustion  chamber.  The  magnitude  of  the  pressure  oscillations  within  the 
combustion  chamber  i3  recorded  along  with  the  phase  angle  difference  between 
the  pressure  and  exit  area  oscillations.  The  chamber  gas  dynamics  are 
predicted  by  an  analytical  model  allowing  the  solution  of  the  remaining 
unknown,  i.e.,  the  propellant  pressure- coupled  response. 

Another  externally  driven  system  employs  a  traditional  acoustic 
measurement  device,  the  impedance  tube.  Impedance  tube  burners  u3e  acoustic 


drives  to  generate  a  standing  wave  in  a  combustion  tube.  The  structure  of  the 
standing  wave  is  measured  by  pressure  transducers  mounted  at  various  points 
along  the  length  of  the  impedance  tube.  The  equation  that  describes  the 
standing  wave  is  derived  from  the  conservation  equations  applied  within  the 
impedance  tube  and  is  a  function  of  the  propellant  admittance.  The  admittance 
i3  deduced  by  incorporating  the  pressure  measurements  into  the  wave  equation. 

All  of  the  aforementioned  methods  result  in  an  indirect  measurement  of 
the  propellant  response  function.  The  pressure-coupled  response  is  derived 
from  a  combustion  analysis  in  conjunction  with  oscillatory  pressure 
measurements.  The  indirect  measurements  are  only  as  accurate  as  the  analyses 
on  which  they  are  based.  A  direct  measurement  of  propellant  pre33ure-coupled 
response  does  not  require  analytical  modeling  to  produce  re3ult3.  One 
technique  utilizing  direct  measurement  is  the  microwave  burner  method.  This 
method  measures  the  transient  burning  rate  of  solid  rocket  propellants  by 
reflecting  a  microwave  signal  from  the  burning  propellant  surface.  The 
oscillatory  burning  rate  can  produce  the  pressure-coupled  response  when 
combined  with  the  simultaneous  oscillatory  pressure  measurement.  A  major 
drawback  to  all  of  the  experimental  techniques  is  that  they  are  unable  to 
produce  accurate  response  function  measurements  at  high  frequencies  (>10  kHz). 
The  problem  of  high  frequency  combustion  instability  has  previously  been 
averted  by  the  use  of  metallic  powders  which  provide  excellent  high  frequency 
damping.  However,  in  the  move  toward  smokeless  propellants,  it  has  become 
necessary  to  remove  the  metal  particulates  and  the  problem  of  propellant 
response  at  high  frequencies  has  reappeared. 

A  magnetic  flowmeter  is  theoretically  capable  of  obtaining  a  direct 
measurement  of  the  ga3  velocity  above  a  burning  solid  propellant  surface  with 
a  frequency  response  greater  than  20  kHz.  The  operation  of  a  magnetic  flow 
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meter  is  based  on  the  principle  first  discovered  by  Faraday  in  which  a 
conducting  gas  (or  any  medium)  moving  through  a  magnetic  field  will  generate 
an  electrical  potential  which  is  normal  to  both  the  magnetic  field  and  the 
flow  velocity.  The  potential  which  can  be  measured  as  a  voltage  will  be 
proportional  to  the  magnitude  of  the  flow  velocity  and  to  the  magnetic  field 
strength.  The  measurement  of  the  oscillatory  gas  velocity  above  a  burning 
solid  propellant  surface  subjected  to  pressure  oscillations  constitutes  a 
direct  measurement  of  the  propellant  pre3sure-coupled  admittance  and  is  the 
basi3  of  the  experimental  portion  of  this  program. 

The  admittance  function  i3  defined  a3  the  nondimensional  ratio  of  the 
combustion  ga3  velocity  to  the  oscillatory  pressure, 

U'/a 

Ab  -  (1) 

P'/P 

where  primes  denote  fluctuating  quantities  and  superscript  bar3  denote  a  mean 
quantity.  The  oscillatory  velocity,  U^’,  is  nondimen3ionalized  by  the  sonic 
velocity  within  the  combustion  chamber,  a.  The  pressure-coupled  response  is 
the  nondimensional  ratio  of  the  ma33  flow  oscillation  to  the  oscillatory 
pressure. 


«„  •  — 

P’/P 

The  response  and  admittance  functions  are  related  by  the  equation, 

Ab  +  Mb  *  T  Mb  Rb 


(2) 


(3) 


Where  Mb  *3  the  Mach  number  of  the  product  gases  as  they  leave  the  propellant 
surface  and  Y  i3  the  ratio  of  specific  heats.  A  magnetic  flowmeter  is  capable 


of  directly  measuring  the  oscillatory  velocity,  Ub' ,  which  in  conjunction  with 
the  oscillatory  pressure,  p',  simultaneously  measured  by  a  pressure  transducer 
produces  the  direct  measurement  of  the  propellant  admittance. 


The  magnetic  flowmeter  operates  according  to  Faraday’s  law  which  states 
that  an  electromotive  force,  emf,  induced  in  a  conductor  moving  relative  to  a 
magnetic  field  is  proportional  to  the  rate  at  which  the  magnetic  field  lines 
are  being  passed  (Ref.  3).  The  average  emf  is  defined  in  equation  form  as: 


emf 


A  <p 
At 


where  A<J>  is  the  change  of  magnetic  flux  relative  to  the  moving  conductor  in  a 
time  interval  At.  In  the  case  of  a  flowmeter,  the  moving  conductor  is  either 
a  gas  or  a  liquid  and  the  change  of  magnetic  flux  is 


A<j>  =  BLUAt  (5) 

where  B  is  the  magnetic  flux  density  (Telsa),  U  is  the  flow  velocity  (m/s), 
and  L  is  the  distance  (m)  the  emf  is  being  measured  across.  Therefore,  the 
induced  emf  expressed  in  volts  becomes 


V  -  UBL  (6) 

A  calibrating  constant  is  required  to  account  for  losses  due  to  nonuniform 
magnetic  fields  and  other  inherent  loss  mechanisms  encountered  in  the  actual 
application  of  the  magnetic  flowmeter. 


V  =  ctUBL  (7) 

The  calibration  factor,  a,  is  nondimensional  and  can  be  determined  either 
theoretically  or  experimentally.  The  operating  principle  of  the  magnetic 
flowmeter  is  presented  schematically  in  Figure  1.  Although  Equation  7  shows 
that  the  potential  is  independent  of  flow  conductivity,  there  exist 
conductivity  limits  for  which  the  equation  holds  true.  A  highly  ionized  flow 
could  allow  magnetohydrodynamic  forces  to  distort  the  velocity  profile.  The 


Princide  of  operation  of  a  magnetic  flowmeter.  Induced 


electric  field  i3  proportional  to  both  flow  velocity  and 


magnetic  field  strength 
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lower  limit  for  conductivity  is  approximately  10~3  mhos/m,  below  which  the 
resistance  of  the  conductor  becomes  as  large  as  the  voltage  measuring  device. 
Experimental  measurements  of  conductivity  in  the  product  gases  of  solid 
propellants  have  been  documented  for  various  propellant  formulations  (Refs.  ^ 
and  5).  All  of  the  formulations  displayed  conductivity  levels  exceeding  the 
magnetic  flowmeter  lower  limit  of  10“3  mhos/m. 

Both  the  mean  and  unsteady  interior  ga3  velocities  of  solid  propellant 
rocket  motors  were  measured  by  Micci  and  Caveny  (Refs.  6  and  7)  using  a 
magnetic  flowmeter.  The  magnetic  flowmeter  produced  well-defined  velocity 
measurements  for  an  AP/HTPB  composite  propellant  in  a  sub3cale  rocket  motor. 
Longitudinal  acoustic  wave3  were  generated  in  the  rocket  motor  by  the 
modulation  of  the  exhaust  nozzle  area.  The  oscillatory  velocity  data 
permitted  the  deduction  of  the  propellant  pressure-coupled  response  from  a 
linear  analysis  of  the  wave  structure.  Holme  (Ref.  8)  investigated  the 
possibility  of  using  magnetic  velocimetry  for  in-flight  measurement  of  rocket 
exhaust  velocity.  The  magnetic  flowmeter  technique  was  tested  in  this  study 
for  several  liquid  propellant  combinations  using  an  alternating  magnetic  field 
to  eliminate  the  d.c.  effects  of  chemical  polarization. 


RESEARCH  PROGRAM 


PRESSURE  COUPLED  RESPONSE  MEASUREMENTS 

The  magnetic  flowmeter  assembly  is  schematically  depicted  in  Figure  2. 

Thi3  figure  presents  the  major  components  necessary  to  measure  surface 
admittances.  Pressure  oscillations  are  produced  by  the  pressure  modulation 
system  (A)  which  i3  located  at  the  converging  nozzle  orifice.  The  pressure 
oscillations  and  the  oscillatory  velocity  are  simultaneously  measured  as  the 
propellant  3trand  burns  by  the  measurement  station  3hown  in  the  cut-away 
section  (B).  The  large  external  magnet  (C)  provides  a  large  d.c.  magnetic 
field  for  the  velocity  measurement  obtained  by  the  voltage  electrodes  at  the 
measurement  station.  The  combustion  chamber  (D)  was  designed  to  consume  small 
amounts  of  propellant  (4.0  cm3)  for  each  run.  The  propellant  samples  were 
formed  into  cylindrical  strands  with  a  diameter  of  1.2  cm  and  a  length  of  3.5 
cm.  The  small  propellant  requirement  makes  the  burner  very  economical  to 
investigate  a  wide  variety  of  propellants. 

Although  commercial  magnetic  flowmeters  use  an  alternating  field  in  order 
to  separate  the  flowmeter  potential  from  chemical  polarization  potentials,  a 
permanent  magnet  was  chosen  for  the  magnetic  flowmeter  burner.  A  permanent 
magnet  is  well  suited  for  the  measurement  of  an  alternating  velocity  field 
because  any  polarization  effects  (which  are  low  frequency)  can  be  filtered  out. 
The  u3e  of  a  permanent  magnet  also  reduced  the  complexity  of  the  burner  by  not 
requiring  additional  circuitry  to  produce  a  magnetic  field.  The  magnet  was  an 
ALNICO  5  permanent  magnet  which  created  a  magnetic  field  with  a  gap  density  of 
approximately  2100  gauss.  The  burner  was  constructed  out  of  304  stainless 
steel  which  permits  the  magnetic  field  to  pass  through  with  very  little 
di stortion. 

The  30nic  nozzle  orifice  was  modulated  by  a  small  toothed  gear  mounted  on 


a  permanent  magnet  d.c.  motor.  The  motor  speed  was  controlled  by  a  Variac 
adjustable  autotransformer  on  an  a.c.  power  line.  The  a.c.  output  of  the 
autotransformer  was  oonverted  to  a  d.c.  voltage  for  operating  the  motor.  This 
arrangement  permitted  steady  performance  at  frequencies  ranging  from  2000  to 
20,000  Hz. 

Figure  3  presents  a  schematic  of  the  data  acquisition  system  utilized  for 
the  magnetic  flowmeter  burner.  The  oscillatory  pressure  signal,  p’,  was 
measured  by  a  PCB  model  105A13  pressure  transducer  with  a  gain  of  100  produced 
by  the  PCB  model  WA  amplifier.  A  high  pass  active  filter  was  used  to  remove 
all  signal  noise  below  1000  Hz.  The  signal  wa3  then  split  for  recording  and 
used  a3  a  reference  signal  by  the  vector  voltmeter  lock-in  amplifier.  The 
lock-in  amplifier,  ITHAC0  model  393*  received  the  velocity  3ignal  directly 
from  the  magnetic  flowmeter  electrodes  and  isolated  the  velocity  component 
oscillating  at  the  same  frequency  as  the  reference  pressure  signal.  The 
output  of  the  lock-in  amplifier  was  the  magnitude  of  the  velocity  oscillation 
in  phase  and  out  of  phase  with  the  pressure  oscillation.  These  two 
measurements  generated  the  real  and  imaginary  parts  of  the  propellant 
admittance  function.  The  recorded  oscillatory  pressure  was  digitized  and 
stored  on  an  IBM  PC. 

A  cold  flow  simulator  was  constructed  to  evaluate  the  magnetic  flowmeter 
burner  design  in  terms  of  the  magnitude  of  pressure  oscillation  and  the 
acoustic  characteristics  of  the  combustion  chamber.  Also,  cold  flow 
simulation  was  necessary  to  determine  any  effects  the  RTV  rubber  diaphragms 
had  on  the  transducer  frequency  response.  The  simulator  duplicated  the 
interior  dimensions  of  the  burner  with  nitrogen  gas  forced  through  a  porous 
metal  plug  to  simulate  flow  from  the  burning  propellant  surface  at  pressure 
levels  up  to  1M.0  MPa.  The  same  motor  and  gear'  pressure  modulation  system  to 
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be  used  on  the  burner  was  applied  in  the  cold  Plow  simulation  to  determine 
whether  It  could  produce  pressure  oscillations  of  sufficient  amplitudes  at 
frequencies  up  to  20  kHz.  The  theoretical  acoustic  characteristics  were 
determined  by  approximating  the  chamber  as  a  close  ended  cylinder. 

The  performance  of  the  pressure  modulation  system  during  cold  flow 
simulation  is  presented  in  Figure  4.  The  pressure  oscillation  amplitudes 
produced  in  the  cold  flow  simulator  are  plotted  as  a  percentage  of  mean 
chamber  pressure,  P,  at  an  excitation  frequency  of  20  kHz.  The  effectiveness 
of  the  exit  area  modulation  declined  initially  as  the  mean  chamber  pressure 
rose  from  1.38  MPa  to  6.9  MPa.  This  behavior  is  illustrated  by  the  decline  in 
the  oscillatory  to  mean  pressure  ratio  in  Figure  4.  The  pressure  ratio 
remained  essentially  constant  from  6.9  MPa  to  13.8  Mpa  at  a  value  of 
approximately  0.7%.  The  amplitude  of  pressure  oscillation  was  dependent  on 
the  mean  chamber  pressure  and  the  gap  between  the  nozzle  exit  plane  and  the 
teeth  of  the  rotating  gear.  The  optimum  gap  distance  was  found  to  be 
approximately  3.0  millimeters.  A  slight  Increase  in  this  distance  resulted  in 
a  large  reduction  in  the  oscillation  amplitude. 

The  cold  flow  simulator  was  also  used  to  test  the  effectiveness  of  RTV 
rubber  diaphragms  of  different  thicknesses.  These  diaphragms  were  necessary 
to  protect  the  pressure  transducer  during  the  combustion  of  the  propellant. 

The  results  of  these  tests  are  plotted  in  Figure  5.  The  diaphragm  thicknesses 
were  1.5  mm  and  3.0  mm  and  they  were  compared  to  an  unprotected  transducer 
over  a  frequency  range  of  4,000  Hz  to  20,000  Hz.  The  mean  chamber  pressure 
was  held  constant  at  3.4  MPa.  The  diaphragms  caused  an  amplitude  loss  which 
was  dependent  on  the  excitation  frequency,  but  the  phase  angle  between  the 
protected  and  the  unprotected  transducers  remained  zero.  Identical  results 
were  obtained  for  the  two  thicknesses  at  every  frequency.  Based  on  these 


experimental  results  and  design  factors,  the  smaller  diaphragm  thickness  was 
selected  for  use  in  the  magnetic  flowmeter  burner.  Figure  5  also  displays  the 
cold  flow  resonant  modes.  The  observed  amplitude  peaks  at  frequencies  less 
than  16,000  Hz  correspond  to  the  predicted  harmonic  resonant  frequencies  for 
longitudinal  mode  oscillation.  The  first  observed  amplitude  peak  occurred 
near  the  second  harmonic  frequency  for  longitudinal  mode  oscillations  with 
additional  peaks  occurring  near  each  consecutive  resonant  frequency  up  to  the 
sixth  longitudinal  harmonic  of  15,720  Hz.  Tangential  mode  oscillation 
occurred  at  frequencies  above  16000  Hz  creating  resonant  modes  which  were 
combinations  of  tangential  and  longitudinal  oscillations.  These  combination 
modes  were  observed  at  two  predicted  resonant  frequencies  of  17,098  Hz  and 
19,885  Hz. 

The  validity  of  the  response  function  measurements  relies  on  the  accurate 
calibration  of  the  magnetic  flowmeter  for  both  velocity  signal  amplitude  and 
phase.  The  amplitude  calibration  consisted  of  determining  the  calibration 
constant,  a,  in  Equation  A.  This  wa3  achieved  by  a  device  designed  to 
simulate  the  velocity  measurement  in  the  magnetic  flowmeter  burner.  The 
conducting  ga3  flow  in  the  burner  was  replaced  by  a  solid  oopper  rod  with  a 


diameter  of  6.0  mm.  The  conductivity  independence  of  the  magnetic  flowmeter 
equation  permitted  this  replacement  as  long  as  the  gas  flow  conductivity  was 
above  the  limit  of  10"3  mhos/m.  An  electromagnetic  shaker  oscillated  the 
oopper  rod  in  the  magnetic  field  produced  by  the  ALNICO  5  permanent  magnet. 
Two  graphite  brushes  picked  up  the  potential  across  the  rod  induced  by  the 
conductor  moving  through  the  magnetic  field.  The  amplitude  of  the  a.c. 
potential  wa3  measured  by  the  lock-in  amplifier  which  used  a  reference  signal 
provided  by  the  shaker  power  supply.  The  velocity  of  the  rod  was  derived  by 
integrating  the  acceleration  measurement  of  a  PCB  model  305A05  accelerometer 
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mounted  on  the  top  of  the  copper  rod.  The  calibration  data  was  recorded  as 
flowmeter  voltage  veraua  oscillatory  velocity  of  the  calibration  device.  This 
data  was  very  consistent  over  a  wide  range  of  velocity  amplitudes.  The  ratio 
between  velocity  and  voltage  wa3  inserted  into  Equation  4  along  with  the 
measured  magnetic  field  strength,  B,  and  the  electrode  separation,  L,  yielding 
a  calibration  constant  of  2.1. 

The  velocity  phase  calibration  was  required  to  correct  the  phase  shift 
induced  by  the  inherent  capacitance  of  the  measuring  circuit  combined  with  the 
high  ga3  resistance.  The  capacitance  induces  an  appreciable  phase  shift  in 
the  velocity  signal  if  the  resistance  of  the  combustion  gas  products  is  large. 
To  quantatively  determine  the  phase  shift  that  occurred  in  the  course  of  the 
experiments,  the  combustion  gas  resistance  had  to  be  measured.  Direct 
measurement  of  the  gas  resistance  using  an  ohmmeter  created  polarization  of 
the  electrodes,  therefore,  an  indirect  method  was  employed  to  determine  the 
gas  resistance.  In  this  method  the  resistance  was  experimentally  determined 
by  measuring  the  phase  shift  induced  on  an  oscillatory  signal  passed  through 
the  combustion  gas  products.  The  advantages  of  this  arrangement  are  that  it 
utilized  the  same  equipment  conf iguration  as  the  magnetic  flowmeter  burner  and 
did  not  require  the  replacement  or  alteration  of  the  existing  combustion 
chamber.  Excellent  run  to  run  agreement  was  achieved  for  the  gas  resistance 
measurements  obtained  at  several  frequencies  of  oscillation,  with  the  mean 
value  of  ga3  resistance  being  40,000  ohm3.  The  value  for  the  gas  resistance 
produced  a  maximum  phase  shift  in  the  velocity  measurement  of  35.6  degrees  at 
a  measurement  frequency  of  20,000  Hz. 

All  propellant  tested  was  of  80S  AP  and  20%  HTPB  composition.  The  first 
response  function  data  wa3  obtained  for  a  20  pm  unimodel  AP  particle  size  at  a 
modulation  frequency  of  8000  Hz.  Figures  6a  and  6b  plot  the  real  and 
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imaginary  parts  respectively  of  the  pressure-coupled  response  for  this 
propellant  formulation.  The  data  is  plotted  as  a  function  of  distance  above 
the  burning  propellant  surface.  The  maximum  distance  plotted,  1.6  mm, 
corresponds  to  the  diameter  of  the  velocity  measurement  electrode.  Mean 
pressures  for  the  te3ts  ranged  from  4.5  to  8  MPa  and  are  listed  for  each  test 
on  Figures  6  through  11.  At  the  pressures  shown  in  Figures  6a  and  6b,  the 
propellant  burn  rate  was  0.58  cm/3.  The  real  part  of  the  response  in  Figure 
6a  displays  excellent  reproducibility.  The  sampling  rate  wa3  increased  to 
obtain  better  spatial  resolution  for  testing  of  a  200  ym  AP  particle  size 
propellant.  Figures  7a  and  7b  show  the  read  and  imaginary  parts  respectively 
of  the  pressure-coupled  response  at  a  modulation  frequency  of  4000  Hz.  The 
value  for  the  read  part  remains  fairly  constant  while  the  imaginary  part 
displays  a  change  in  value  through  the  flame  zone.  Figures  8a  and  8b  3how  the 
real  and  imaginary  parts  at  a  modulation  frequency  of  8000  Hz.  For  these 
tests  the  sampling  rate  was  increased  to  its  then  maximum  obtainable  value  of 
250  Hz.  At  thi3  frequency,  the  real  part  Of  the  response  increases  with 
distance  above  the  propellant  surface  while  the  imaginary  part  remains 
constant  at  a  slight  negative  value.  At  the  chamber  pressures  shown  in 
Figures  7  and  8  the  200  ym  AP  propellant  burn  rate  was  approximately  0.7  cm/s. 

Additional  pressure- coupled  measurements  were  made  by  two  undergraduate 
students  for  their  senior  thesis  project.  Measurements  were  made  at 
modulation  frequencies  of  12,  16,  and  20  kHz  for  the  200  ym  unimodal  AP 
propellant.  Difficulties  were  encountered  obtaining  a  sufficiently  high 
si gnal-to-noise  ratio  and  the  data  exhibited  greater  scatter  than  at  the  lower 
modulation  frequencies.  Figures  9a  and  9b  plot  the  real  and  imaginary  part3 
respectively  of  the  pressure-coupled  response  at  a  modulation  frequency  of 
12,000  Hz.  The  data  for  the  two  run3  at  thi3  frequency  show  agreement  for  the 
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Fjg.  8a  Real  part  of  pressure-coupled  response  as  a  function  of 

distance  above  propellant  3urface  for  200  um  AP  propellant 
tested  at  8000  Hz. 


(1.0  L 


2.0 


l«siw»v  P*n  of  Ru»o#sr  Fvkctioh 
8000  Hj 
700*"  »P 


•  TO4S824  (1000  PSIA) 
c  TOS0827  (3000  PSIA) 


•••••• 

lo 

°  o 


-2.0  lb  0 


-4.0  b 


oCooo°®(,*?*»pcc  - 
0°  °  • 


oo o®°0o°_  -50*** 


S®cr 


•  •••  nuu 

■'o  o  *•  0°oo*o  #0# 

Or,  C  ••••  0 


I _ I _ I _ I _ I _ I _ 1 _ I - 1 

o  0.4  O.E  1.2  1.6 

Dis'«»a  A*ovt  P»pi>tu»NT  Si'*f*ct  («*) 

fig*  8b  Imaginary  part  of  pressure-coupled  r'e3ponse  as  a  function  of 
distance  above  propellant  surface  for  200  ym  AP  propellant 


tested  at  8000  Hz. 


Real  Pa»t  or  Rumdse  Function 

12.000  Hi 

20Q*  AP 


O  PC720206  (1000  PSIA) 
•  PC 7502 06  (1000  PSIA) 


10  \  Og-  *08 


•  0*0  o  o 

•  O  _ 


Rg  o  - -8&-  -  - - c "Ofo*'***" "  •*“ 

o  ••  o 

•  o° 

-10  *  *  oo  0o 


Fig.  9a 


DISTANCE  ABOVE  PtOFILLANT  SURFACE  0*0 

Real  part  of  pressure-coupled  response  as  a  function  of 
distance  above  propellant  surface  for  200  ym  AP  propellant 
tested  at  12.000  Hz. 
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Fig.  9b  Imaginary  part  of  pre33ure-coupled  response  a3  a  function  of 
distance  above  propellant  3urface  for  200  um  AP  propellant 
tested  at  12,000  Hz. 


real  part,  however,  agreement  cannot  be  claimed  for  the  imaginary  part. 

Figures  10a  and  10b  show  the  real  and  imaginary  parts  respectively  of  the 
pre33ure-coupled  response  at  a  modulation  frequency  of  16,000  Hz.  Agreement 
between  the  two  runs  made  at  this  frequency  was  not  obtained  for  either  the 
real  or  imaginary  parts.  Figure  11  plots  both  the  real  and  imaginary  parts  of 
the  pressure-coupled  response  for  the  one  run  made  at  a  modulation  frequency 
of  20,000  Hz.  Depletion  of  the  supply  of  200  pm  unimodal  AP  propellant 
prevented  additional  testing  at  these  high  frequencies.  At  the  chamber- 
pressures  shown  in  Figures  9  through  11,  the  propellant  burn  rate  was 
approximately  0.7  cm/3. 


PRESSURE-COUPLED  RESPONSE  ANALYSIS 

The  effort  to  obtain  an  analytical  model  for  predicting  the  high 
frequency  solid  propellant  response  functions  follows  the  method  U3ed  by  Tien 
(Ref.  9).  This  method  numerically  integrates  a  fifth  order  system  of  ordinary 
differential  equations.  The  flame  is  assumed  to  be  one-dimen3ional ,  premixed, 
laminar,  and  a  one-3tep  forward  chemical  reaction  is  taking  place.  The 
specific  heats,  molecular-  weights,  and  coefficient  of  heat  conduction  are 
taken  to  be  constant,  and  the  Lewis  nianber  is  equal  to  unity.  The  governing 
equations  for  the  reacting  flow  above  the  30lid  propellant  surface  are 
continuity,  energy,  and  specie  conservation  given  by 
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where  w  is  the  reaction  rate,  q  is  the  heat  of  reaction,  vj  is  a  function  of 
the  reaction  stoichrometric  coefficients  and  the  molecular  weight  of  the 
species,  and  Y  is  the  mass  fraction.  A  perturbation  analysis  was  utilized 
with  a  small  amplitude  sinusoidal  time  oscillation  assumed  for  the  variable 
and  with  both  mean  and  perturbation  quantities  being  functions  of  height  above 
the  burning  propellant  surface.  Since  the  pressure  wave  length  is  always  much 
larger  than  the  flame  thickness,  the  pressure  was  assumed  to  be  a  function  of 
time  only  and  a  momentum  equation  is  not  required.  Tien's  model  i3  the  only 
one  to  assume  an  unsteady  gas  phase  reaction  in  addition  to  the  unsteady  heat 
transfer  to  the  burning  30lid  propellant.  The  boundary  condition  equations 
resulted  from  solving  the  energy  equation  in  the  solid  phase  and  by  assuming 
an  isentropic  flow  very  far  from  the  flame. 

A  fourth  order  Runge-Kutta  method  was  utilized  to  integrate  from  the 
outer  flame  edge,  x  -  •,  to  the  propellant  surface,  x  -  0.  The  integration 
was  carried  out  in  thi3  direction  to  allow  a  linear  conversion  method  to  be 
used  to  calculate  the  boundary  values.  The  linear  conversion  method  separated 
the  solution  into  homogeneous  solutions  and  a  particular  solution  (Ref.  10). 
Each  of  the  homogeneous  solutions  was  expressed  in  tern>3  of  the  x  -  •  boundary 
value.  The  dependent  variable  for  each  homogeneous  solution  was  given  an 
initial  boundary  value  of  unity  and  a  value  of  zero  for  the  particular 
solution.  The  system  of  equations  was  then  integrated  for  each  set  of 
boundary  values.  The  resulting  solutions  were  summed  and  placed  in  the 
boundary  equations  to  determine  the  oorrect  boundary  values  of  x  «  •. 

Figure  12  plots  the  mean  temperature  and  reaction  rate  as  a  function  of 
distance  above  the  burning  propellant  surface  which  is  at  location  x  -  0. 

These  curves  were  obtained  by  integrating  the  mean  energy  equation  and  are 
required  as  input  for  the  integration  of  the  perturbation  equations.  The 


results  given  by  the  unsteady  gas  phase  model  or  T'ien  were  then  compared  with 
the  pressure-ooupled  admittances  measured  for  the  200  pm  unimodel  AP  composite 
propellant  at  4000  and  8000  Hz.  The  model  was  utilized  to  calculate  the 
oscillating  velocity  component,  u',  versus  height  above  the  burning  propellant 
surface  through  the  flame  zone  for  comparison  with  experimentally  measured 
values.  The  two  key  parameters  which  were  varied  are  the  nondimensionalized 
frequency,  to,  and  the  ratio  of  gas  residence  time  to  the  30lid  residence  time, 
Z.  Figure  13  shows  a  representati  ve  result  for  w-2  and  E-0.0122.  It  can  be 
seen  that  the  model  predicts  a  real  component  of  the  fluctuating  velocity  that 
levels  out  at  a  constant  value  with  increasing  distance  from  the  propellant 
surface  while  the  imaginary  component  steadily  becomes  more  negative  with 
increasing  distance.  This  compares  with  the  experimental  data  3hown  in  Figure 
7.  Z  and  u>  were  then  varied  to  determine  if  the  T'ien  model  could  reproduce 
the  functional  dependence  of  the  admittance  on  the  frequency  of  the  pressure 
oscillations.  Figure  14  shows  the  real  part  of  the  admittance  for  two 
frequencies  a  factor  of  two  apart.  It  can  be  seen  that  for  the  Z  of  0.024*1 
the  higher  frequency  produces  a  higher  real  part  of  the  admittance.  This  was 
also  observed  experimentally  as  can  be  seen  by  comparing  Figures  7  and  8. 
However,  the  model  still  predicted  an  imaginary  part  with  a  constant  negative 
slope  unlike  the  measured  value  for  the  imaginary  part  shown  in  Figure  8  at  a 
frequency  of  8000  Hz.  It  is  felt  at  this  time  that  the  T'ien  model  is  too 
simplistic  in  the  chemical  reaction  kinetics  to  accurately  model  all  the 
observed  behavior  of  the  solid  propellant  admittance. 

VELOCITY-COUPLED  RESPONSE  MEASUREMENTS 

A  burner  to  measure  solid  propellant  velocity-coupled  responses  was 
designed  based  on  the  knowledge  gained  from  operation  of  the  pressure-coupled 


idmittance  aa  a  function  of  nondimenslonal  distance 
surface  for  two  values  of  the  frequency,  u. 


burner.  Here  the  velocity-coupled  response  is  defined  to  be  the  response  of 
the  propellant  burning  rate  to  an  oscillatory  cross-flow  velocity  adjacent  to 
the  propellant  surface.  Figure  15  shows  the  propellant  geometry  and  the 
orientation  of  the  electrodes  and  magnetic  field  to  permit  the  simultaneous 
measurement  of  the  cros3-flow  oscillating  velocity,  V',  and  the  oscillatory 
velocity  of  the  ga3  leaving  the  propellant  surface,  Ub'.  It  is  a  2-D  3lab 
motor  with  two  pairs  of  velocity  measuring  electrodes  embedded  in  two  facing 
propellant  3lab3.  Three  of  the  electrodes  penetrate  one  propellant  slab  with 
the  fourth  electrode  in  the  opposing  slab.  One  pair  of  electrodes  will 
measure  the  cross-flow  velocity  above  the  propellant  surface  (V')  and  the 
other  pair  will  measure  the  velocity  of  the  hot  combustion  ga3se3  leaving  the 
propellant  surface  (Ub').  Correlating  the  two  measurements  will  give  the 
velocity-coupled  response.  Figure  16  shows  the  velocity- coupled  magnetic 
flowmeter  burner.  A  is  the  burner  assembly,  B  is  the  permanent  magnet,  C  is  a 
cutaway  view  of  the  propellant  slabs  and  velocity  electrodes  and  D  is  the 
rotating  toothed  gear  used  to  modulate  the  burner  orifice  and  generate  a 
standing  wave  in  the  combustion  chamber.  The  velocity  measurement  station  is 
located  at  the  midpoint  of  the  chamber  where  a  velocity  antinode  and  a 
pressure  node  would  be  situated  for  a  standing  wave  at  or  near  the  chamber 
fundamental  longitudinal  frequency.  This  was  done  in  order  to  maximize  the 
oscillatory  oro33-flow  velocity  to  obtain  a  high  signal-to-noi3e  ratio  and  to 
minimize  pressure  oscillations  which  would  generate  a  pressure-coupled 
response.  The  electrodes  are  rigidly  mounted  and  the  height  of  the 
measurement  station  above  the  burning  propellant  surface  increases  as  the 
propellant  regresses,  ranging  during  a  burn  from  0  to  9  mm  at  burnout.  Figure 
17  is  a  detailed  drawing  of  the  velocity-coupled  magnetic  flowmeter  burner 
showing  the  locations  of  the  mean  chamber  pressure  transducer,  the  electric 
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15  Slab  propellant  geometry  and  orientation  of  electrodes  and  magnetic 
field  for  the  almultaneoua  measurement  of  the  oscillatory  cro33-flow 
velocity,  an<f  We  oaclll atory  velocity  of  the  gas  leaving  the 
propellant  surface,  ly. 


Fig.  16  Veloci ty -coup led  magnetic  flowmeter  burner  uaed  to  make  3imultaneou 
measurement  of  oscillatory  cro33-flow  velocity,  V1,  and  oscillatory 
velocity  of  gaa  normal  to  propellant  surface,  Uh*.  Components: 
A-burner  assembly,  B -permanent  magnet,  C -cutaway  view  of  propellant 
slabs  and  velocity  electrodes,  D-rotating  toothed  gear . 


match  used  for  ignition  and  the  over pressurization  vent  at  the  head  end  of  the 
chamber  and  the  sonic  nozzle  at  the  opposite  end.  Velocity-coupled  response 
data  is  currently  being  obtained  with  this  device. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


There  are  two  major  shortcomings  for  the  methods  presently  used  to 
measure  solid  propellant  pre33ure-coupled  response  functions.  They  are  unable 
to  produce  reliable  results  at  high  frequencies  (>  10  kHz)  and  the  majority  of 
the  methods  rely  on  an  indirect  means  of  response  function  measurement.  One 
of  the  objectives  of  this  study  was  to  develop  an  apparatus  that  utilizes  a 
magnetic  flowmeter  to  obtain  direct  measurements  of  the  pressure-coupled 
response  functions  for  solid  propellants.  The  device  was  required  to  perform 
these  measurements  over  a  frequency  range  up  to  20,000  Hz  while  consuming 
small  amount  of  propellant.  A  burner  apparatus  was  designed  and  constructed 
to  meet  these  specifications.  A  magnetic  flowmeter  was  incorporated  in  the 
apparatus  in  order  to  measure  the  velocity  fluctuations  of  the  combustion  gas 
products  occurring  during  unsteady  combustion  of  solid  propellant  samples. 
This  velocity  measurement  was  the  essential  feature  of  the  magnetic  flowmeter 
burner  because  it  permits  the  pressure-coupled  response  function  to  be 
obtained  directly  from  the  experimental  data. 

The  magnetic  flowmeter  burner  provided  pressure-coupled  response 
measurements  for  two  formulations  of  AP/HTPB  composite  propellant  at  selected 
frequencies.  The  AP  particle  size  differed  between  the  two  formulations,  20 
pm  AP  and  200  pm  AP.  The  magnetic  flowmeter  provided  oscillatory  velocity 
measurements  as  the  electrodes  were  exposed  to  the  combustion  ga3  products. 
Mean  and  oscillatory  pressure  measurements  were  recorded  simultaneously  with 
the  oscillatory  velocity  measurements.  These  measurements  produced  the 
surface  admittance  function  for  the  burning  propellant  samples.  A  simple 
relationship  produced  the  pressure-ooupled  response  function  from  the 
admittance  function.  Repeatable  behavior  was  noted  in  the  reeil  and  imaginary 
parts  of  the  pressure-coupled  response  for  the  20  pm  AP  and  200  pm  AP 


propellant  samples. 

The  validity  of  the  response  function  measurements  relies  on  the  accurate 
calibration  of  the  magnetic  flowmeter.  The  magnitude  of  the  oscillatory 
velocity  signal  wa3  readily  calibrated  by  the  data  generated  by  the  velocity 
calibration  device.  This  data  wa3  consistent  over  a  wide  range  of  velocity 
amplitudes.  The  velocity  phase  shift  correction  wa3  found  to  be  dependent  on 
the  combustion  gas  resistance.  The  resistance  was  experimentally  determined 
by  measuring  the  phase  shift  induced  on  an  oscillatory  signal  passed  through 
the  combustion  gas  products.  The  advantages  of  this  arrangement  are  that  it 
utilizes  the  same  equipment  configuration  a3  the  magnetic  flowmeter  burner  and 
did  not  require  the  replacement  or  alteration  of  the  existing  combustion 
chamber.  Excellent  run  to  run  agreement  was  achieved  for  the  gas  resistance 
measurements  obtained  at  several  frequencies  of  oscillation. 

The  magnetic  flowmeter  burner  satisfied  the  objective  of  this 
experimental  study  by  producing  direct  pressure-coupled  response  function 
measurements.  The  burner  proved  capable  of  measuring  pressure-coupled 
response  at  frequencies  beyond  the  limits  imposed  on  other  methods.  The  small 
amount  of  propellant  consumed  by  the  magnetic  flowmeter  burner  for  each  run 
makes  the  technique  very  economical  for  extensive  testing  of  various  solid 
propellants.  Some  difficulty  was  experienced  with  obtaining  the  measurement 
of  the  oscillatory  velocity  at  a  frequency  other  than  the  resonant  frequency. 
This  problem  did  not  seriously  hamper  the  effectiveness  of  the  magnetic 
flowmeter  burner  and  can  most  likely  be  remedied  by  minor  design 
modifications. 

The  external  method  of  inducing  pressure  oscillations  within  the 
combustion  chamber  of  the  burner  produced  smaller  than  expected  amplitudes  and 
the  oscillatory  pressure  signal  was  subjected  to  d.c.  shifts  due  to 


temperature  effects  and  partial  exhaust  nozzle  blockage.  A  high  pass  filter 
and  an  100  gain  amplifier  was  incorporated  into  the  circuit  to  condition  the 
oscillatory  pressure  signal.  The  improved  pressure  signal  was  inadequate  for 
the  lock-in  reference  3ignal.  The  small  pressure  oscillations  resulted  in 
velocity  measurements  with  a  low  signal-to-noise  ratio.  Measurements  at  an 
excitation  frequency  of  8000  Hz  were  of  sufficient  amplitude  to  be  readily 
distinguished  from  background  noise.  This  was  due  to  the  excitation  frequency 
being  near  the  first  longitudinal  resonant  frequency  of  the  combustion  chamber. 
Cold  flow  simulation  has  shown  that  the  oscillation  amplitude  is  much  larger 
at  the  chamber  resonant  frequencies.  However,  the  velocity  signal  was  of  the 
same  order  of  magnitude  a3  the  noise  at  frequencies  other  than  resonant 
frequencies.  The  velocity  signal  was  harder  to  discern  for  the  burner  tests 
at  other  frequencies.  The  test  had  to  be  repeated  several  times  to  verify  the 
data.  A  method  to  increase  the  velocity  signal  or  reduce  the  amount  of 
background  noise  is  needed  before  effective  measurements  can  be  taken  at 
other  frequencies. 

Magnetic  flowmeter  burner  effectiveness  can  be  improved  by  modifying  the 
chamber  design  to  allow  a  larger  diameter  propellant  strand.  Additional 
propellant  consumption  would  be  offset  by  the  benefits.  A  larger  exit  nozzle 
orifice  size  would  accomodate  the  larger  diameter  propellant  strands  tested 
over  the  same  range  of  mean  chamber  pressures.  Nozzle  blockages  by  combustion 
residue  will  decrease  and  mass  flow  through  the  nozzle  will  increase.  The 
greater  mass  flow  permits  higher  amplitude  pressure  oscillations  generated  by 
the  pressure  modulation  system.  In  addition  to  an  improved  reference  signal 
for  the  lock-in  amplifier,  the  larger  oscillations  within  the  combustion 
chamber  will  increase  the  velocity  measurements  by  the  magnetic  flowmeter. 

Another  benefit  of  the  larger  diameter  propellant  strands  involves  the 
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effect  of  the  velocity  electrodes  embedded  in  the  strands.  The  propellant 
burning  surface  deforms  around  the  electrodes  as  they  are  exposed  by  the 
receding  propellant.  The  calibration  of  the  flowmeter  velocity  measurements 
a33umes  an  uniform  flow  vector  from  a  planer  burning  3urfaee.  Thi3  assumption 
13  rendered  invalid  if  large  deformities  alter  the  average  propellant  3urface 
level.  The  effect  of  the  deformities  caused  by  the  electrodes  would  be  le33 
significant  on  a  larger  burning  3urface  area. 

The  unsteady  gas  phase  analysis  of  T'ien  has  been  used  to  successfully 
predict  the  trends  of  the  pre33ure-coupled  response  a3  a  function  of  frequency. 
The  model  numerically  integrates  the  governing  equations,  continuity,  energy 
and  specie  conservation,  to  predict  the  mean  and  oscillating  components  of  the 
flow  variables:  velocity,  temperature,  and  specie  mass  fractions.  Numerical 
results  and  experimentally  measured  pressure-coupled  response  values  are 
compared.  However  a  more  accurate  combustion  process  model  i3  necessary  to 
quantitatively  predict  the  propellant  high  frequency  pressure-coupled 
response. 

Techniques  validated  in  the  pressure-coupled  magnetic  flowmeter  burner 
were  applied  to  design  and  construct  a  high  frequency  velocity-coupled 
magnetic  flowmeter  burner.  Two  pairs  of  velocity  electrodes  simultaneously 
measure  oscillatory  cros3-flow  velocity  and  oscillatory  ga3  flow  velocity 
normal  to  the  propellant  surface.  The  geometry  utilized  was  a  two-dimensional 
3lab  motor.  The  30nic  nozzle  is  modulated  to  generate  a  fundamental  mode 
standing  wave  in  the  chamber  with  the  measurements  station  located  at  the 
velocity  anti  node.  Experimental  data  i3  currently  being  taken  to  determine 
the  ability  of  the  magnetic  flowmeter  technique  to  measure  a  velocity-coupled 


response. 
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